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Generalization of the Euler Angles 


Malcolm D. Shuster 1 and F. Landis Markley 2 


Abstract 

It is shown that the Euler angles can be generalized to axes other than members of an 
orthonormal triad. As first shown by Davenport, the three generalized Euler axes, hereafter: 
Davenport axes , must still satisfy the constraint that the first two and the last two axes be 
mutually perpendicular if these axes are to define a universal set of attitude parameters. 
Expressions are given which relate the generalized Euler angles, hereafter: Davenport angles , 
to the 3-1-3 Euler angles of an associated direction-cosine matrix. The computation of the 
Davenport angles from the attitude matrix and their kinematic equation arc presented. The 
present work offers a more direct development of the Davenport angles than Davenport’s 
original publication and offers additional results. 

Introduction 

The Euler angles [ 1-4] are defined as the angles of a sequence of three rotations 

R(h v n', n"; *>, rf, VO = R(A 3 , #)R( n l5 <p ) , (1) 

which represent the attitude, in particular, the attitude of a rigid body. Here, n 1# 
n 2 , and are selected from the set {1, 2, 3}, where 



V 


"o' 


V 

1 = 

0 

, 2 = 

1 

, and 3 = 

0 


0 


0 


1 


In general, we denote column vectors by bold sans serif letters. A caret here 
denotes a unit column vector. The primes denote that the column vectors are each 
representations with respect to a different abstract basis (as seen by an inertial 
observer). In this case, the basis is the current basis of the body-fixed coordinate 
system, which changes (from an inertial point of view) as the body rotates. 
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In order that the representation in terms of the Euler angles have the required 
three degrees of freedom, we must further stipulate that 

n x ^ n 2 and n 2 ^ n 3 (3) 

Given this restriction, there are twelve possible sets of Euler angles: six symmetric 
sets, whose labels are written as 

1-2-1 1-3-1 2-3-2 2-1-2 3-1-3 3-2-3 

and six asymmetric sets, designated by 

1-2-3 1-3-2 2-3-1 2-1-3 3-1-2 3-2-1 

In each label the first (leftmost) integer denotes the first rotation axis. For example, 
the 1-3-2 set of Euler angles correspond to n x = 1, n 2 = 3, and n 3 = 2. The 
asymmetric sets have been called variously Cardan angles, Bryant angles and Tait 
angles. 

Regarded as the representation of body-fixed axes, the Euler axes are generally 
chosen from a right-hand orthonormal triad of column vectors. That this orthonor- 
mal triad need not be limited to the representations of the body coordinate axes 
with respect to themselves (that is, to the set {1,2,3}) should be obvious but 
will be demonstrated rigorously below. What we wish to know in particular is 
whether it is possible to construct a representation of the attitude in terms of Euler 
rotations about three arbitrary non-orthogonal axes as seen from the body-frame. 

In the present note we prove that a universal representation of the attitude in 
terms of three consecutive Euler rotations about arbitrary non-orthogonal axes is 
not possible. However, we shall show that an extension of the definition of the 
Euler angles does indeed exist. The Euler angles for those sets of Euler axes 
which accommodate the representation of any attitude we shall refer to as universal 
attitude parameters , or as a universal representation , in order to distinguish them 
from the angles about three axes for which some attitudes cannot be represented. 

Universality of the Conventional Euler Angles 


To the best of our knowledge, all texts present formulas for extracting the 
conventional Euler angles from the rotation matrix, but none demonstrate rigorously 
that the Euler angles can represent an arbitrary rotation matrix. The proof is quite 
simple, and we offer it here. 

We particularize our discussion to the 3-1-3 set of Euler angles. We shall show 
later in this report that any of the remaining eleven sets of conventional Euler 
angles can be obtained from the formula for extracting the 3-1-3 Euler angles. 
Hence, it the first step will be to prove that the 3-1-3 Euler angles are a universal 
parameterization of the attitude. 

The explicit evaluation of equation (1) for n : = 3, n 2 — 1, and n 3 = 3 yields 3 

R 313 (^ d, VO EE *(3, 1P)R( i, d)R(3, ip) 


cip ap — siped s<p 
—sipap — cip cd sip 
sd sip 


cip sp + sip cd c<p sip sd 
—sip s(p + cip cd ap cip s'd , 
—sd ap cd 


( 4 ) 


3 Henceforth, we shall discard the primes on the Euler axis representations with the understanding that 
these are always with respect to current body axes. 
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where op = cosy?, sip = siny>, etc., and we have written the axis indices as 
subscripts on R. The elements of the attitude matrix (direction-cosine matrix) are 
the components of the initial body axes with respect to the final body axes. In order 
to demonstrate that the parameterization of R 313 can realize any proper orthogonal 
matrix, it is sufficient to show that the representation of the coordinate axes of the 
initial coordinate system with respect to themselves (namely, 1, 2 and 3) can be 
transformed into any other right-hand orthonormal triad of column vectors, i.e., the 
representation with respect to an arbitrary set of arbitrary right-hand orthonormal 
coordinate axes. 

Thus, we examine first 


©3 = ^313 (^» i 3 ) 3 — 


sin d sin ip 
sin d cos ip 
cos'd 


( 5 ) 


Clearly, d and (n/2-ip) are the spherical angles of an arbitrary unit vector. Hence, 
by a suitable choice of d and ip, e 3 can be made to coincide with an arbitrary unit 
column vector, in particular, the third member of the target right-hand orthonormal 
triad. We will assume that d and ip have been so chosen. 

Examine next 


©1 = ^313 1 


cos ip cos ip — cos d sin ip sin ip 
— sin ip cos ip — cos d sin ip cos ip 
sin d sin ip 


( 6 ) 


This last column vector can be written as 

e x = cosipu l + sinypu 2 , (7) 



cos ip 


— cos sin ip 

Ul = 

— sin ip 

. “2 = 

— cos d cos ip 


0 


sin d 


Obviously, u x and u 2 are each unit column vectors and 

Ui-e 3 = U2-e 3 = u 1 -u 2 = 0. (9) 

The vectors Oj and u 2 span the plane perpendicular to e 3 , and ip can be chosen 
to generate any unit vector in that plane. Since the first member of the target 
right-hand orthonormal triad must lie in that plane, there is a value of <p for which 
e 1 coincides with that unit column vector. Thus, for suitable choices of the three 
Euler angles, R 313 will transform 3 and 1 into any pair of mutually perpendicular 
unit column vectors. 

Finally, because both the initial column vectors { 1, 2, 3}) and the target column 
vectors {e 1 ? e 2 , e 3 } are each right-hand orthonormal sets and R 313 (ip, d, ip) is 
proper orthogonal by construction, it follows that 


^313 (<p, VO 2 — ^313 (v>, VO (3 x i) 

= (^'313 (ip, ■&, vo 3) x (r 313 (ip, ip) i) 

~ © 3 X ©1 = ©2 


(10a) 

(10b) 

(lOcd) 
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which completes the proof. Therefore, the 3-1-3 Euler angles are a universal 
representation of the attitude. The proof can be repeated with minor modifications 
for any of the other eleven conventional sets of Euler angles, but the universality 
of these remaining sets will follow as an corollary of the result for the generalized 
Euler angles. 

Generalized Euler Angles (Davenport Angles) 

Consider a three-parameter rotation given by 


R{h v n 2 , n 3 ; ip, i>) = R( n 3 , V>) R( n 2 , #) R(n v <p) , (11) 

where now, there are no restrictions on n x , n 2 , and n 3 , except that 

Hi 7 ^ ^2 ^2 ^3 ' ( 12 ) 

Can any rotation be represented in this way? A necessary condition for this to be 
true is that any unit column vector can be transformed into any other unit column 
vector by this sequence of rotations. For the generalized Euler angles to be a 
universal representation of the attitude, this condition must be satisfied. 

Consider now the transformation 


v = R(n 1 , n 2 , n 3 ; <p, ip)^. (13) 

Then 

n 3 • v = n 3 R(n 2 , «?) n x . (14) 

We shall show that equation (14) cannot be satisfied for arbitrary n 2 , n 2 , n 3 , and 

v. 

Writing Euler’s formula as [3] 


R(n, 0 = I 3x3 + sin C [[ n ]] + (1 -cos<) [[ n ] 


(15) 


with 


0 


IUJJ S I -u, 0 u, | , (16) 

L u 2 — «! 0 J 

it follows straightforwardly that the condition on d becomes 

n 3 ■ v = (n 3 • n,) + n 3 • (n 2 x (n 2 x nj) - siml (n 3 ■ (n 2 x n,)) 

- cosd (n 3 ■ (n 2 x (n 2 x n^)) . (17) 


Define now 

B cos ot = — n 3 • (n 2 x (n 2 x n 1 )) and B sina = -n 3 • (n 2 x n 1 ) , (18) 

or, equivalently, 

a = arctan 2 [h x ■ (n 2 x n 3 ), -n 2 ■ (n 2 x (n 2 x fi 3 ))] , (19a) 

> 0 , (19b) 


B = \J [n 3 • (n 2 x fij)] 2 + [n 3 • (n 2 x (ii 2 x n x ))] 2 



Generalization of the Euler Angles 


5 


where arctan 2 (^ x ) is the function which yields the arc tangent of y/x in the 
correct quadrant. This corresponds to the function ATAN2 in the FORTRAN 
programming language. In terms of these new variables, equation (17) becomes 

n 3 • v = 0 + B cos($ - a ) , (20) 

where 

0 = (n 3 • fii) + n 3 • (n 2 x (n 2 x nj) = (n 3 • n 2 )(n 2 ■ nj . (21) 

The right member of equation (20) can assume any value between 0 - B and 
0 + B. Therefore, a solution will exist for if and only if 


0 + B>n 3 -v>0-B. (22) 


However, since n 3 • v can assume any value between -1 and + 1, it follows from 
equation (22) that 0 and B must satisfy 


B> 1-0 and £>! + /?. 

(23) 

Thus, we require that 

£> l + |0|. 

(24) 

On the other hand, defining 

u = n 2 x n x , 

(25) 

it follows that 


B 2 = [n 3 • u ] 2 + [n 3 • (n 2 x u )] 2 


= |u | 2 n 3 { u u T + (n 2 x u)(n 2 x u) T } n 3 . 

(26) 

Now, fi 2 and u are orthogonal. Hence, n 2 , u, and n 2 x u, form 
triad, and therefore 

an orthonormal 

n 2 n^ + u u T + (n 2 x u)(n 2 x u) T = / 3x3 . 

(27) 


It follows that 

B 2 = |u | 2 rtj [/ 3x3 - n 2 n^ ] n 3 = |u | 2 (1 - (n 3 ■ n 2 ) 2 ] = |n 3 x n 2 | 2 |n 2 x Ail 2 . (28) 
We have thus 


/? = (n 3 • n 2 )(n 2 • fii) , B = |n 3 x n 2 | |n 2 x fijl . (29ab) 

Equations (24) and (29) can be satisfied simultaneously if and only if 

B = l and 0 = 0. (30ab) 


Hence, we require that 

Hj JL n 2 and n 2 ~L ng (3 lab) 

as a necessary condition that the generalized Euler angles be able to represent an 
arbitrary attitude. Equations (31) were first discovered by Paul Davenport [5 ], who 
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also proved sufficiency, as we shall below. For this reason we will refer henceforth 
to the generalized Euler angles as the Davenport angles. 

Equation (31) is a less restrictive condition than the generalization of equa- 
tions (2) and (3) to an arbitrary right-hand orthonormal triad, which would have 
required further that (n 1 ■ n 3 ) be either 0 or 1. The additional degree of freedom 
allowed by equation (31) is the angle between A : and n 3 . 

We shall now prove that the condition expressed by equation (31) is sufficient 
for the set of Davenport axes to permit a universal parameterization of the attitude 
in terms of Davenport angles. Let us write 

n 3 = R{h 2i A) , (32) 

where h l and n 2 are orthogonal. Then, 

n 3 = cosAA 1 -sinA(A 2 x A x ), -7r<A<7r, (33) 

which clearly satisfies n 3 1 n 2 and is the most general column vector satisfying 

this condition. A is the angle from A x to n 3 , defined to be positive in the 

counter-clockwise direction about n 2 . The symmetric sequences of Euler angles 
correspond to A = 0, while the asymmetric sequences correspond to either A = n/2 
or A = — 7r/2. If n 2 is the cyclic follower of n ls then A = n/2. If it is the anticyclic 
follower, then A = -7r/2. 

It follows [3] that 

R( n 3 . #) R(n v <p) 

= R(R( n 2 , A)n 1( V>) #(n 2 , tf) R(n u <p) (34a) 

= R( n 2 , X)R{n 1 , i>)R T { n 2 , A).R(n 2 , , 0)R(n l , <p) (34b) 

= R(n 2 , \)R{h v n 2) ft,; <p, tf', V), (34c) 

where 

- A . (35) 

If A is an arbitrary proper orthogonal matrix, then wc wish to find angles (</?, ip) 
which satisfy 

A = R( fij, n 2 , n 3 ; y>, VO ~ R(n 2 , A) R{ n v n 2 , A x ; <p, t?', VO . (36) 

To show that this is possible, let C be the proper orthogonal matrix which satisfies 

Cit 1 = 3 , and C n 2 = 1 ■ (37) 

Since h 1 and n 2 are orthogonal, the matrix C exists and is given by 

C = [ n 2 (A x x A 2 ) A x ] , (38) 

where the expression for C is given as the transpose of a proper orthogonal matrix 
labeled by its column vectors. It follows from equations (36) and (38) that the 
triplet (v?, if, VO must satisfy 

R( 3, 1, 3; ip, f) = R t (1, A )CAC t . 


(39) 
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The triplet (</?, tf', ip) is now simply the 3-1-3 set of Euler angles representing the 
matrix which is the right member of equation (39). Since the 3—1—3 Euler angles 
have been shown to be a universal representation of the attitude, we know that a 
solution always exists. This proves the sufficiency of equations (31). Q.E.D. 

As a by-product of our proof we have that 

R(n v n 2 , n 3 ; <p, i>) = C T R{\, A) R(3, 1 , 3; <p, 0- A, $)C (40) 

for any set of axes satisfying equation (31) and C given by equation (38). As an 
immediate corollary we have that any of the twelve conventional sets of Euler angles 
is a universal attitude representation. For a suitable choice of A and C, expressions 
similar to equation (40) can be obtained for all twelve sets of conventional Euler 
angles. 

Recall that A is a function solely of the axes. 

A = arctan 2 ((n 1 x n 2 ) ■ n 3 , n x ■ n 3 ) . (41) 

This formula also assumes that equation (31) holds. 

We remark that while the three Davenport axes must satisfy equation (31) in 
order for the representation to be universal, the parameterization of the attitude 
can still be useful when only equation (12) is satisfied, if one knows a priori that 
the axes permit a realization of the attitude matrix for the values of interest. 

Extracting the Davenport Angles 


To determine the Davenport angles from a given direction-cosine matrix we note 
the relationships 


njAhx = cos($ - A) , 

(42a) 

n 2 A hj = sin(i9 - A) sin ip , 

(42b) 

(n 2 x n 3 ) T A hi = - sin(t9 - A) cos ip , 

(42c) 

n 3 A n 2 = sin(r? - A) sin <p , 

(42d) 

n^A (n x X n 2 ) = — sin(t9 - A) cos <p . 

(42e) 

Except for the fact that has been replaced by - A, the right members of 
equations (42) are each identical within a sign to the elements of the direction- 


cosine matrix given in equation (4). 

From equation (42a) we have immediately 

t3 = A + arccos(n^A nj , (43) 

If we choose the principal value of the arc cosine, then will be single valued and 
lie in the range 

A < d < A + 7T . (44) 

For A < < A + 7r, so that sin(tf - A) > 0 the two remaining Davenport angles are 

given by 


tp = arctan 2 [n 3 . A h 2 , — n 3 A (n 1 x n 2 )] , 
ip = arctan 2 [n^4 fij , -(n 2 x n 3 ) T /4n 1 ] . 


(45a) 

(45b) 
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The cases where sin($ - A) = 0 are treated in the next section. 

Singularity of the Davenport Angles 

The representation in terms of Davenport angles must become singular when 


^1 — ^^3 » ( 46 ) 

which follows from an examination of the resulting equations 
#(n 3 , ip) R(h 2 , 0) R(n v <p) 

= R( n 3 , ip) R( n 2 , t9) <p) R T ( n 2 , t9) R(h 2 , tf) (47a) 

= -R(A 3 , rp) R(± n 3 , tp) R(h 2 , tf) (47b) 

= R(h 3 , i>±<p) R{ n 2 , t?) (47c) 

= R(n 2 , i9) R(n l , p±ip) , (47d) 


where we assumed equation (46) in going from equation (47a) to equation (47b). 
It is easy to see from equation (42a) that a singularity in the Davenport angles 
occurs when sin (# - A) = 0. 

Equations (45) are inappropriate for the computation of p and ip at a singularity, 
when d - A = 0 or 7r, because all of the arguments of the arctan 2 functions will 
vanish. In that case we must turn to four other “elements” of the direction-cosine 
matrix, namely, 


A T * ^ 

n 2 A (n x x n 2 ) = sin tp cos ip + cos tp sin ip cos(i9 — A) = a , 

(n 2 x n 3 ) T A n 2 = cos tp sin ip + sin tp cos ip cos($ - A) = b , 

n 2 A n 2 = cos tp cos ip - sin tp sin ip cos(tf — A) = c , 

(n 2 x n 3 )^A (fli x h 2 ) = sin tp sin ip — cos tp cos ip cos(tf — A) = d . 

From these relationships it follows that 

[I ± cos($ — A)] sin(<£ ± ip) = a±b, 

[1 ± cos(i9 — A)] cos (tp ± ip) = c d , 

from which it follows that 


(48a) 

(48b) 

(48c) 

(48d) 

(49a) 

(49b) 


p + ip — arctan 2 (a + 6, c — d) for cos($ - A) ^ -1 , (50a) 

tp — ip = arctan 2 (a - 6, c + d) for cos(tf — A) ^ +1 . (50b) 

The.se equations provide a more accurate means for calculating either tp + ip or 
p — ip (but, unfortunately, not both) when the attitude is very close to a singularity 
than do equations (45). 
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The Davenport angles extracted by this procedure and that of the previous section 
will not always agree with the conventional Euler angles when a conventional Euler- 
axis set is used. This is due to trivial angular ambiguities of 27r and to the two-fold 
ambiguity in the Davenport angles, which is related to the well-known two-fold 
ambiguity of the conventional Euler angles. The general relation for the two-fold 
ambiguity in the Davenport angles can be obtained from equation (40) and the 
similar relation for the 3—1—3 Euler angles, namely, 

R( 3, i , 3; <p, d, Xl>) = R{3, i, 3; p + tt , -tf , (51) 

with the result 

n 2 , n 3 ; <p, ip) = R(n v n 2 , n 3 ; <p + n, 2A - ■d, ip - ir) , (52) 

which summarizes the result for the 3-1-3 Euler angles (A = 0) and the 3-1-2 
Euler angles (A = 7r/2). The specification of the formulas of this section and the 
previous section to the extraction of the conventional Euler angles is left as an 
exercise for the reader. 

Kinematics 


The kinematic equation for the Davenport angles is identical to the general 
expressions for the conventional Euler angles. 4 We may write [3] 

« = V>n 3 +tf#(n 3 , ip) n 2 + <pR{n 3 , i>)R{ n 2 , (53) 


where is the body-referenced angular velocity vector. This may recast in the 
form 

<?] \<P 

= R(n 3 , xp)S(n 1} n 2 , n 3 ; t9) = M(<p, ip) , (54) 

ip ip 

with ^(n^ n 2 , n 3 ; tf) represented in terms of column vectors as 
S(n l5 n 2 , n 3 ; tf) = [R{n 2 , i9)n x | n 2 | n 3 ] . 


(55) 


Simple forms exist for the matrices n 2 , n 3 ; tf) and M(<p y d, ip) as explicit 

functions of the conventional Euler angles [4]. Those for the Davenport angles 
are more complicated. 

More useful for simulation is the inverse of M(ip, ip), which satisfies 


d 


= M 1 (c^, d, ip) oj , 


whence 

■d, ip) = S _1 (ni, n 2 , n 3 ; d) R T ( n 3 , ip ) , 

with 


S' _1 (n 1 , n 2 , n 3 ; ■d) = 


1 

(#(n 2 , «?)n x ) • (n 2 x n 3 ) 


[n 2 x n 3 ] T 
[n 3 x R(h 2 , i?)n!] T 
.[(«(n 2 , tJ)n,) x n 2 ] 7 ’. 


(56) 

(57) 


( 58 ) 


l The kinematic equation for the Davenport angles was also derived by Davenport [5]. 
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The right member of equation (58) can be reduced to 


D 


l n 2 x n 3 J 

[cost? (n 3 x n x ) - sin'd(h l x n 2 ) x n 3 ] 
[costfi^ x n 2 - sini9 n l ] T 


where the denominator D is given by 

D = (R( n 2 , ) n L ) • (n 2 x n 3 ) = - sin(tf - A) 


(59) 


(60) 


Discussion 

We have given a presentation of the Davenport angles in the same detail that 
is normally given to the standard twelve sets of Euler angles. The connection has 
been made between the Davenport angles and the 3-1-3 set of Euler angles, and 
this has been used to develop further relationships of the Davenport angles. A 
corollary of our work is that if a reorientation of a spacecraft is to be accomplished 
by rotational maneuvers about three axes, then these must be Davenport axes. An 
application of the Davenport angles to attitude estimation is given in Reference 6. 
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